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Abstract 

A catalogue of all non-isomorphic simple connected regular matroids M. of cardinality n < 15 
is provided on the net. These matroids are given as binary matrix matroids and are sieved 
from the large pool of all non-isomorphic binary matrix matroids of cardinality < 15. For 
each A4 its Tutte polynomial is determined by an algorithm based on internal and external 
base activity. 
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We assume familiarity with basic matroid concepts such as contraction matroids, dual matroids, 
or the flat lattice [B]. A famous binary matroid is the Fano matroid Fj on the set {a, b,- ■ ■ ,g} 
which can be defined as the column matroid of this 3x7 matrix over GF{2): 



A 



It is essential in the sequel that two binary matroids are isomorphic if and only if some ob- 
viously sufficient condition holds for any two respective matrix representations. Specifically, for 
any column representation of a binary n-matroid, every other column representation is obtained 
by replacing the rows rj of the matrix by a set of equivalent rows, i.e. generating the same 
subspace of GF(2)"', and by subsequently permuting the labelled columns [6, 10.1.3]. Thus, in 
our case we could take this representing matrix for Ff. 
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The fact that this, like every matrix representation of F^, has the same column set as A is par- 
ticular to Ff. It makes it easy to decide whether or not a 3 x 7-matrix over GF{2) represents F^. 
It is well known that the dual M (here = Ff) of a representable matroid A4 is obtained by 
taking the orthogonal complement Y = X-^ oi X = rowspace(yl) with respect to the canonical 
scalar product r • s in GF{2) , and letting A' be any matrix whose rows constitute a basis of Y . 
For instance, one verifies that ri ■ Sj = Q for all 1 < z < 3 and 1 < j < 4: 
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Therefore, the column matroid of A' must be isomorphic to F^. In fact, whether or not 
any 4x7 matrix of rank 4, with columns labelled a,b, . . . ,g from left to right, yields a column 



matroid isomorphic to Fj, is most easily checked by testing whether its rows are orthogonal to 

The simplification 7W of a matroid AA is the matroid formed by removing all loops and 
identifying any multiple points into a single point. Recall that a matroid is regular if it is 
representable over every field. There are various ways to characterize the regular matroids among 
the class of binary matroids. For us the algorithmic most convenient one was the following. 



Theorem 1 [6, 13.4. l(ii)]: A binary matroid M of rank k is regular if and only if there 

id 
7 • 



neither is a rank k — 3 fiat U with Ai/U ~ Fj, nor a rank k — A flat V with Ai/V ~ Ff 



In view of Theorem 1 the following procedure suggests itself to generate up to isomorphism 
all regular matroids up to size N: 

a) Generate up to isomorphism all binary matroids Ai up to size A^. 

b) For each fixed A4 of rank k < N compute all rank k — 3 flats Ui,- ■ ■ ,Us and all rank A; — 4 
flats Vi,---,Vt. 



c) Check whether M./Ui ~ Fj for some 1 < i < s, or M./Vj ~ Fj for some 1 < j < f . 

Here comes the section break up. The details of a), b), c) are discussed in sections 2, 3, 
4 respectively. Section 5 outlines a novel algorithm for computing the Tutte polynomial of a 
matroid. It was used to calculate the Tutte-polynomials for all regular matroids generated. 

2 Generating binary matroids up to isomorphism 

Expanding upon the introductory remarks, each binary matrix matroid Ai of cardinality n and 
rank k can be represented by a binary k x n-matrix A of rank k. This representation, however, 
is not unique, since all matrices of the form G ■ A for G € GLfc(2), the group of all invertible 
k X k matrices over GF{2), represent the same matroid. 

Two matrix matroids represented by the k x n-matrices A and B over GF{2) of rank k are 
isomorphic if and only if there exists an n x n-permutation matrix P and some G € GLfc(2) so 
that B = G ■ A- P. Thus, the isomorphism class of the matrix matroid given by A is the orbit of 
A under the action of GL,fc(2) x Sn on the set of all k x n-matrices of rank k. This orbit consists 
of all matrices of the form G ■ A ■ P where G € GLfe(2) and P is a permutation matrix. 

Next we consider these matrices A as functions from n := {1, . . . , n} to GF{2) , where A{iy , 
the transposed of A{i), is the i-th column of the matrix A. The multiplication A- P oi matrices 
is then replaced by composing A with a suitable permutation vr G 5,1. Thus we have a group 
action of GLfc(2) x Sn on the set of all functions from n to GF{2)^ determining matrices of rank 
k. Since the symmetric group Sn acts on the domain of these functions, it is not important in 
which position a vector v G GF{2)^ appears in A. We are only interested in how often v appears 
in A. The group action of GLfc(2) x S„ can be replaced by an action of GLfc(2) on the set of 
all jS^-orbits of functions mapping n to GF{2) . The S'n-orbit of a function (or matrix) A is a 
multiset of elements of GF{2)^ of size n. It contains each column of the matrix A together with 
its multiplicity. Hence, for obtaining the isomorphism classes of binary matroids it is enough 
to consider the action of GLfc(2) on the set of all multisets of elements of GF{2)^ of size n. 
Such multisets are written as functions / : GF{2)^ -^ Z>o := {0,1,2,...} where f{v) is the 
multiplicity of the column v G GF{2)^ in /. Therefore, we call / the multiplicity function of a 
multiset of GF{2)'^ of size n. Consequently, the multiplicity functions / of a multiset of GF{2)^ 
of rank k and size n satisfy the following two properties: 

a) There exist k linearly independent vectors v G GF(2) such that f{v) > and 



b) E^,GGF{2)* /(^^ 



n. 



Finally, we introduce the function p : GF{2)^ -^ {0, 1, . . . , 2*^ — 1} for labelling the elements 

V = (fi, . . . ,Vk) S GF{2)'^ with nonnegative integers. (Thus p is a rank function on GF{2)^ .) 

The label of v is 

k 

p{vi,...,Vk) :=^Vj2^~^, 

3=1 

where we identify the elements and 1 of GF[2) with the corresponding integers. The unit 
vectors (1, 0, . . . , 0), (0, 1, 0, . . . , 0), . . . , (0, . . . , 0, 1) in GF{2f are then labelled by the powers 
of 2, i.e. by 1, 2, ... , 2^~^ , respectively. 

Conversely, given a label r E {0,1,..., 2 — 1} we determine the unrank function p^^ : 
{0, 1, . . . , 2^^ — 1} — )• GF{2)^ by p^^(r) = v = {vi, . . . , Vk) which is obtained from the binary 
representation of r as r = Ylj=i ^i^-'^^ with Vj G {0, 1}. 

Using p, we determine a permutation representation GL,fc(2) of GLfc(2) on the set {0, 1, ... , 
2^ - 1}. Consider B G GLfc(2) and j G {0, 1, . . . , 2'^ - 1}, then ttb defined by iTBiJ) := 
p{B ■ (p"^(j))^) is a permutation of {0, 1, ... ,2'' - 1}. 

Using p we consider the multiplicity functions / of multisets of GF{2) of size n as functions 
/ from {0, 1, . . . , 2^^ — 1} to Z>o with the property X^j=o /O) ~ "-• ^^ write these functions 
/ as vectors / = (/(O), /(I), . . . , f{2^ — 1)). The set of these vectors is totally ordered by the 
lexicographic order. The group GLfc(2) acts in a natural way on the domain of these functions. 
This induces a group action on the set of all functions from {0, 1, . . . , 2^^ — 1} to Z>o. We consider 
the lexicographic largest element of the GLfc(2)-orbit of / as the standard representative of this 
orbits. 

The functions / are supposed to describe matrix matroids of rank fc, whence, these matrices 
contain a set of k linearly independent vectors. Due to the labelling of the vectors of GF{2)^ 
and since the standard representative is the largest vector in its orbit, a standard representative 
contains the labels of the k unit vectors in GF{2)^ with nonzero multiplicities. Thus for i G 
{2-' I < j < A;} we have f{i) > 0. Moreover, it is easy to see that standard representatives / 
satisfy f{i) > f{r) for all r > i and all i G {2-' | < j < k}. 

In order to determine a complete list of representatives of binary matrix matroids of size n 
and rank k we apply Read's method of orderly generation (cf. [3 H [21 [5] ) in the following setting. 
The group GLfc(2) acts on the domain of all multiplicity functions / from {0, 1, . . . , 2^ — 1} to 
Z>o which satisfy X^»=o fU) — ''^- We are listing all these functions lexicographically decreasing 
starting with the lexicographically largest one. Due to our construction we restrict ourselves 
to those / with f(i) > and f{i) > f{r) for all r > i and i G {2^ | < j < k}. Each of 
these functions is tested whether it is the lexicographically largest element in its GLfc(2)-orbit. 
If so, it is considered to be a canonic representative of a binary matrix matroid. Otherwise it 
represents a matrix matroid which is isomorphic to a matrix matroid which was already found 
before. In this case it is sometimes possible to obtain information which multiplicity function 
/ should be tested in the next step (for more details see [3j). In several situations we do not 
proceed with the lexicographically next / but we are allowed to do some jumps in the set of all 
these functions. 

A matrix representation A = Af of a multiplicity function / associated with a matrix 
matroid of rank k and size n is then obtained from /, by building a, k x n- matrix consisting 
of /(i) columns of {p^^{i)) for i G {0,1,..., 2 — 1}. It is easier to write down just the 
labels of the columns of this matrix, thus we obtain a vector rj := (ri, . . . ,r„) of labels with 
fi, < Ti^+i for i^ < n and {{ly \ I < ly < n, r,^ = i}\ = f{i) for i G {0, 1, . . . , 2^^ — 1}. Due to 
this construction, if /i and /2 are multiplicity functions of size n and /i is larger than /2 with 



respect to the lexicographic order, then the vector rj^ is smaller than rj^ in the lexicographic 
order of all vectors in {0, 1, . . . , 2^^ — 1}". Since we present the binary matroids in form of these 
vectors (ri,...,r„), the standard representative of an isomorphism class is now the smallest 
element in its orbit and the representatives are listed monotonically increasing according to the 
lexicographic order. 

All these computations were done in SYMMETRICA [8]. For obvious reasons we restrict 
ourselves to loopless matroids. Thus the columns of the matrix matroids belong to GF{2)^ \ {0}, 
therefore, /(O) = for all /. 

This way we obtain complete lists of representatives of the isomorphism classes of all loopless 
binary matrix matroids and of all simple loopless binary matrix matroids of size < 15 and rank 
< 7. 

For example, there are two simple loopless binary matrix matroids of rank 3 and size 4. They 
are given by 

■ 1 1 
1^1,2,3,4) which represents the matrix I 1 1 



and by 



1 

10 1 

(1,2,4,7) which represents the matrix | 1 1 

11 



These matroids were found from the list of all binary matroids of rank 3 and size 4. It con- 
tains the three vectors ri := (1,1,2,4), r2 := (1,2,3,4), r^ := (1,2,4,7). The unit vectors 
have labels 1, 2, 4. Since the rank k = 3, the multiplicity functions are functions from 
{0,1,..., 7} to Z>o satisfying X^j=o/0) — ^- The corresponding multiplicity functions are 
/i = (0, 2, 1,0, 1, 0, 0, 0), /2 = (0, 1, 1, 1, 1, 0, 0, 0), h = (0, 1, 1,0, 1, 0, 0, 1) and the correspond- 
ing multisets of GF{2f of rank 3 and size 4 are Ai = {(1, 0,0), (1,0, 0), (0, 1, 0), (0, 0, 1)}, 
^2 = {(1,0,0), (0,1,0), (0,1,1), (0,0,1)}, and ^3 = {(1,0,0), (0,1,0), (0,0,1), (1,1,1)}. 

3 Getting the flats of rank A: — 3 and k — A 

Suppose 7W is isomorphic to the column matroid of the /c x n matrix A over GF{2). The circuits 
of M are well known to be the inclusion-minimal members among the sets supp(2;) := {z € n | 
Zi = 1} where z ranges over the subspace (rowspace(A))-'-. Actually, circuits come up only in 
section 5. For the time being we rather need the (say) h many cocircuits Di of M (i.e. the 
circuits of Ad ) which can be found as the minimal members among the sets supp(z) where 
z ranges over the 2'^-element rowspace of A. As for any matroid (representable or not), the 
complements Hi := n\Di {1 < i < h) are precisely the hyperplanes of Ai, i.e. the rank k — 1 
flats. 

The rank k — 2 flats T clearly are the inclusionwise maximal sets among the sets HiOHj (1 < 
i < j ^ h). Similarly the rank A; — 3 flats U are obtained from the T's, and then the rank /c — 4 
flats V in the same manner from the C/'s. 

4 Calculations of contractions 



Let A4 be a binary matroid and U some flat such that the simplification A4/U of the contraction 
Ai/U is isomorphic to Fj. If / C [/ is any generating set, i.e. its closure cl{I) is U, then also 
Ai/I = Fj. It turns out that independent generating sets I '^U are the most appropriate since 
matrix representations of Ai smoothly transform into matrix representations of Ai/I. 



Example: Each binary matroid Ai on the set {a, b,- ■ ■ ,m, n}, with basis B = {a, b, c, d, e}, can 
be represented by a matrix of type 



(^5,^) 



where Is is the 5x5 identity matrix. Provided the rows are labelled by the base elements in 
the right order, this matrix can without loss of information be replaced by the standard matrix 
A of Ai with respect to the basis B: 
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For instance, one still reads off from A that f = a + b + e. It turns out that when the independent 
set / happens to be a subset of B, then B — I is a basis of Ai/I, and the standard matrix A* of 
A4/I with respect to the basis B — I is obtained from A by cancelling the rows with labels in /. 
But what if the independent generating set /' of the flat U is not a subset of B? Then we 
extend /' to some other basis B' of A4. To fix ideas, say /' = {b,g} ^ B. We switch g with any 
element in B that features in the column of g, say with c. By pivoting with respect to Xc,g = 1 
the matrix A transforms to some matrix A': 
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Generally, the matrix C" = {nj^s) arising from C = (x^y^s) by pivoting with respect to Xa,i3 / 
is defined [Ox, p. 84] by 



2/7,<5 := 



a^7,<5 + x^^i3Xa,s, if 7 7^ a and 6^/3; 
x.y^s, otherwise. 



It turns out that A' above is the standard matrix of Ai with respect to the basis B' := {B — 
{g}) U {c} = {a,b,g,d,e} of Ai. Since now /' C B', the standard matrix A* of Ai/I' with 
respect to the basis B' — I' = {a, d, e} is 
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From this one reads off tliat h, i are parallel in M./I' and that M-jV ~ F-j. Therefore, M. is not 
regular by Theorem 1. D 

We mention that it pays off to handle at once all rank A: — 3 flats \J oi M that "fit" the 
current standard matrix A oi M with respect to base B. Namely, they are the (,) many flats 
U = cl{I), where / ranges over all {k — 3)-element subsets of B. 

5 Computing the Tutte polynomial 

For an element a of a matroid Ai we denote hy A4 — a and A4/a the deletion and contraction 
of a respectively. Usually the Tutte polynomial T{M.,x,y) of Ad is calculated via the deletion- 
contraction formula 



(1) 



T{M,x,y) 



T{M-a,x,y) + T{M/a,x,y) 



which holds for all a € Al which neither are isthmuses nor loops. A sophisticated application 
of (1) is discussed in [4J but it only concerns graphic matroids. 

The second author of the present article wrote a very different Mathematica program to 
calculate T{Ai,x,y) for any matroid A4 whose circuits and cocircuits are known. It is based on 
the formula 



(2) 



T{M,x,y) 
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X 



i{B)ye(B) 



from [9, p. 236], where the sum is over all bases i? of A^ and where i(-B), e.{B) are the internal 
and external activities of B. 

Let us provide the gist of the algorithm; a detailed discussion will be given elsewhere. First 
the external activity of B is defined [9^ p. 234] as the number e{B) of x a n\B for which there 
is a circuit C such that 

C Q B U {x} and x = max(C). 

Dually, the internal activity of B is the number i{B) of x a B for which there is a cocircuit D 
such that for B' := n\B we have 

D C B' U {x} where x = max(L'). 

Of course, getting all circuits and cocircuits of A4 is easier said than done. For graphic A4 
a method is given in [lOj . When A4 is provided as the column matroid of a matrix A over a 
finite field K, we may proceed as in section 3. For n < 15 and K = GF{2) this straightforward 
approach suffices. 

As to getting the bases consider the rank 5 polygon matroid Ad on the edge set 8 of the 
graph: 




One reads off that the circuits Cj are (using shorthand notation) 

126, 678, 1287, 3456, 12345, 34578. 

Processing them one by one with the n-algorithm of [11] yields a compact encoding of the 
family Mod of all X C 8 that contain no set Ci. Thus, in our situation. Mod is the family of all 
independent sets of M.: 
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Specifically Mod is given as a disjoint union ri U r2 U rs where each row rj consists of 
certain 0, 1-vectors X of length 8, corresponding to subsets of 8 in the usual way. The string 
of symbols (say) n2n2n2 in ri by definition means that only vectors X are allowed that have 
at least one in a position occupied by an n2- Similarly for uiui and n^n^. Thus ri contains 
(2^ - 1)(23 - 1)(22 - 1) = 63 vectors X, one of them say 01101100. There is a systematic 
(recursive) procedure to determine the number of sets of fixed cardinality in such a multivalued 
row. But here we proceed ad hoc. Namely, the sought bases of being the maximal independent 
sets, each base occurs as a maximal member of some row rj. Here, incidentally, each row- 
maximal member also is a base (since it has cardinality 5). The row-maximal members of rj 
are obtained by filling each n-bubble to full capacity with I's in all possible ways. For instance, 
ninininini in r2 can be filled with four I's in 5 ways. It follows that r2 features 5 • 2 = 10 
bases. Altogether there are 

12 + 10 + 6 = 28 

bases B. It is plausible that for each base B the parameters e{B) required in (2) can now be 
calculated fast. The numbers i{B) can be obtained the same way; after all they are just the 
numbers e{B') for the dual bases B' . 

We mention that the bottleneck in the sketched method is getting the circuits and cocircuits 
(the bases and cobases are then readily found as sketched). How the circuits can be retrieved 
faster than by browsing the whole cycle space, and how the cocircuits arise from the circuits 
(even the chordless ones), is work in progress. 



6 Lists of non-isomorphic regular matroids 

Complete lists with tables of numbers of non-isomorphic regular matroids can be found on the 
web: 



http : //www . uni-graz . at/~f ripert/html/matroids/matroidejieu . html| 



In addition to these tables we also provide complete lists of representatives of these matroids, 
and their Tutte polynomials. 

As it was previously explained, first we determine lists of non-isomorphic binary matroids. 
For 1 < n < 15 and 1 < A; < 7 there are complete lists of the representatives of the isomorphism 
classes of all (a) loopless, (b) simple loopless, (c) connected simple loopless, and (d) connected 
loopless binary matroids of rank k and size n. As the standard representative of a matroid we 
choose the lexicographic smallest representative in its isomorphism class. 

By testing these matroids for regularity, we obtain complete lists of standard representatives 
of the isomorphism classes of (a) connected simple loopless and (b) connected loopless regular 
matroids of rank k and size n again for 1 < n < 15 and 1 < /c < 7. 



The lists of matroids of rank > 8 were computed by changing from a matroid to its dual. If A4 
is a binary matroid of size n and rank k, then its dual, 7W^ is a binary matroid of size n and rank 
n — k. Two binary matroids are isomorphic if and only if their duals are isomorphic. Moreover, 
a matroid is regular if and only if its dual is regular. Using this method we determined complete 
lists of (in general non-standard) representatives of the isomorphism classes of (a) connected 
simple loopless and (b) connected loopless regular matroids of rank A; > 8 and size n < 15. 

Additionally, the Tutte polynomials of all connected simple loopless regular matroids are 
determined and can be downloaded in form of Mathematica Notebooks. 
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